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Abstract 

A method is developed which allows to determine the first-order and the 
second-order magnetoelastic coefficients of a magnetic bulk material from the 
ab-initio calculation of the magnetocrystalline anisotropy energy as function of 
a prestrain €q. Explicit results are given for bcc Fe, and they agree well with 
experimental data obtained from the magnetostrictive stress measurements 
for epitaxial Fe films. 

PACS: 75.70.-i, 75.80.-hq, 71.15.-m 



In recent years magnetic devices based on magnetic films technologies have attracted 
a considerable interest, e.g., magnetooptical recording media or magnetoresistive devices 
based on the giant magnetoresistive and the tunnel magnetoresistive effect designed for sen- 
sors or magnetostrictive random access memories. Thereby the magnetic anisotropy plays 
an important role, for instance, the issue of perpendicular anisotropy for the magnetoopti- 
cal recording or the demand for soft magnetic layers with weak anisotropy as part of the 
magnetoresistive devices. It has been shown by numerous investigations that the magnetic 
anisotropics of magnetic films grown epitaxially on a substrate may strongly deviate from 
those of the respective bulk materials. The reason for this deviation is in general ascribed to 
several different effects. First, there are contributions to the anisotropy originating from the 
free surface of the magnetic layer and from the interface between layer and substrate, as well 
as from the morphology of the film due to a heterogeneous film growth. The influence of all 
these effects must decrease with increasing thickness of the film. What remains for a film of 
thickness larger than typically 10 nm is the effect of the magnetoelastic coupling to the film 
strain induced by the lattice mismatch between film and substrate. Because the epitaxial 
film strain may be of the order of several % which is much larger than typical magnetostric- 
tive strains of 10^^ to 10"'' and because of the dependence of the magnetoelastic coupling 
energy on the direction of the magnetization (see below), this may result in new magnetic 
anisotropics different from that of the unstrained bulk. The numerous experiments on the 
effect of epitaxial strain on the magnetic properties of magnetic films are reviewed in Ref. 
||l[|. An ab-initio study of this effect within the framework of density functional theory is the 
purpose of the present letter. 

To be more specific we consider a material which is cubic in the unstrained state (Fe, 
for instance). Then the density of the magnetoelastic coupling energy may be written (up 
to the second order in the strain e^j, omitting the terms including the shear strain e^, i ^ j, 
which are not required for the situation discussed below) as 
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Here Bi and Dn, D12 represent magnetoelastic coupling coefficients of the first and the 
second order, and the denote the direction cosines of the magnetization referred to the 
cubic axes. For Fe grown epitaxially on a cubic (100) surface which is a prototype system 
experimentally investigated intensively there are epitaxial strains en = €22 which 

may be different for various atomic layers of the film In the following we adopt a 

simple model where we consider only the average strain of the film which depends on the 
thickness of the film, i.e., we assume en = €22 = eo- The strain £33 then may be obtained 
by minimizing the total energy density / = /me + /ei, with the elastic energy density (again 
omitting the terms containing the shear strains): 

/el = ^C*!! {e\^ + + €33) + C12 (eiie22 + £22633 + e33eii) (2) 
and with the cubic elastic stiffness constants Cj,-. This yields 
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_ 2Ci2eo + 5ia3 
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Because the magnetostrictive contribution to the strain originating from /me is much smaller 
than eo we find €33 ~ — 2C12/C11. 

From (1), (3) we obtain the strain dependent part, /mca(eo)5 of the magnetocrystalline 
anisotropy energy density: 

/mca(eo) = /me(eO, «! = !)- /me(eO, "3 = 1) = ^0 + ^iCq + ^269, (4) 

with 

k,= i?i(l + 2Ci2^^), (6) 

fc2 = iZ^n(l-4§^). (7) 

The term is negligible and arises from the fact that we fix the strains en = 622 = 
eo independently on the direction of the magnetization. Furthermore, when changing the 



direction of the magnetization there is a change in the magnetostrictive stress ti = 
according to 

An = ri(ai = 1) - ri(a2 = 1) = + Dueo- (8) 

An experimental determination of Ati (exploiting the change of the bending moment that 
is created by the film onto the substrate ) as function of the layer thickness and hence 

as function of eo then yields the two magnetoelastic coupling coefficients Bi and Du. For 
Fe on MgO (100) Koch et al. H obtained Bi = -3.2 MJ/m^, Dn = 1.1 GJ/m^ (±10%), 
and for Fe on W(IOO) Enders et al. g found Bi = -3 MJ/m^, Du = 1 GJ/m^. The 
values for Bi extracted from the film experiments agree rather well with the bulk value 
of Bi = —3.44 MJ/m^. There are no values of Du obtained from bulk measurements for 
comparison and therefore the two experiments were considered as the first determination of 
the second-order magnetoelastic coupling coefficient of bulk Fe by a film experiment. 

It was pointed out |]1],0,^ that due to the large strains accessible by epitaxial film growth 
the effective first-order coefficient -Befr defined as B^s = Bi+DutQ changes sign from negative 
to positive for Fe with decreasing film thickness, i.e., increasing eo, and this clearly demon- 
strates that the magnetic anisotropy energy depends dramatically on the film thickness, a 
result which is most relevant for the design of the magnetic film devices (see introduction). 
The linear strain dependence of B^s failed to describe the experimental data for film thick- 
nesses below 10 nm, most probably because then the effects of the surface, interface and film 
morphology become relevant. 

In the present paper we determine the magnetoelastic coupling coefficients Bi and Du 
for cubic Fe by the ab-initio density functional theory. To do this, we calculate the mag- 
netocrystalline anisotropy energy density /mca(eo) as function of the strain en = €22 = eo 
imposed to the bulk material, represent the data by a quadratic polynomial in eo according 
to eq. (4) and determine Bi and Du from eqs. (6,7), inserting the elastic stiffness constants 
C12 and Cu which we have also obtained ab initio. 

We have performed the calculations using the WIEN97 |5| code which adopts the 
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full-potential linearized augmented plane- wave (FLAPW) method For the exchange- 
correlation potential the local-spin-density (LSDA) functional by Perdew and Wang and 
the generalized-gradient-approximation (GGA) functional by Perdew et al. |p were used. 
The total energy minimizations on the non-strained bcc Fe gave us the equilibrium lattice 
parameters a = 5.2 ao for LSDA and a = 5.34 ao for GGA where ao denotes Bohr's radius. 
The calculated ratio 2C12/C11 is 1.08 for LSDA and L13 for GGA. The experimental values 
are a = 5.42 ao and 2C12/C11 = 1.17. 

Numerically the most difficult step is the calculation of /mca which is due to the spin-orbit 
coupling (SOC). First, we calculate the self-consistent electronic structure in the scalar- 
relativistic approximation using A^^ k vectors with A^k = 21 in the total Brillouin zone 
(BZ) which correspond to the 762 k vectors in the irreducible part of the Brillouin zone 
(IBZ). The criterion for the self-consistency is the difference in the charge densities after 
the last two iterations being less than 2 x 10~^e/(a.u.)^. The contribution of the SOC is 
determined perturbatively using the second variational method p!0|JTl|] . The quantity /mca 



is calculated by applying the force theorem []T^|I3| as the difference between the sums of 



the perturbed eigenvalues for the different magnetization directions. Fig. 1 represents the 
convergency test for the calculation of /mca with respect to A^k- The data are for the case 



with a = 5.4 ao and c = 5.2 ao, using GGA. The modified tetrahedron [|I4[ and the Gaussian 



smearing pSiJTq] integration schemes were used. The proper convergency with the Gaussian 
smearing was achieved by setting the smearing parameter as F/Ak. The suitable values of F 
for the particular case are roughly from the interval between 6.8eV and lO.leV which follows 
from the curves in Fig. 1. All the final calculations of /mca were performed with A^k = 51 
(17576 k vectors in the IBZ) using both the tetrahedron and the Gaussian smearing method 
in order to minimize the numerical uncertainties. The estimated accuracy is ±l/ieV/unit cell 
marked by the horizontal lines in Fig. 1. 

Fig. 2 shows the calculated magnetocrystalline anisotropy energy density /mca with re- 
spect to the lateral strain eo. The numerical data are well fitted by quadratic polynomials 
as predicted by eq. (4). From the calculated parameters ki and k2 the magnetoelastic 



coefficients Bi and Dn are determined according to eqs. (6,7). Table I summarizes the 
theoretical results in comparison with the experimental data from There is a big dis- 

crepancy between the LSDA result for Bi and the experimental result, whereas the GGA 
result is much closer to the experiment. This is in line with the calculation of the magne- 
toelastic coefficient Aioo of unstrained bulk Fe by Wu et al. who also obtained a strong 
deviation from the experiment when using LSDA but a satisfactory agreement when using 
GGA. The calculated second-order magnetoelastic coupling coefficient Dn for bulk bcc Fe 
matches the experimental value obtained from the measurements on epitaxial thin films very 
well, especially the value from the GGA calculation. The agreement represents the direct 
proof that the experimental results can be really ascribed to the pure strain effect on 

the magnetoelastic properties and that the measurements of the magnetostrictive film stress 
as function of the film thickness can provide the second-order coupling constant Dn of the 
bulk which is hard to obtain by bulk measurements. 

We close with an important warning. The experiments and the present theory 

demonstrate that the magnetoelastic properties of thin epitaxial films may deviate signifi- 
cantly from that described by the ffist-order magnetoelastic coupling coefficients of the bulk 
and that one has to take into account the second-order terms of the bulk. The change of 
the magnetostrictive stress Ari obtained when switching the magnetization from [010] to 
[100] as function of the epitaxial strain eo then may be expressed by an effective ffist-order 
coefficient B^q = Bi + DntQ. However, the correct result for the coefficients /cq, ki and /c2 
of the polynomial expansion (4) of /mca niay ™t be obtained by neglecting the second-order 
terms in eq. (1) and instead replacing Bi in the ffist-order term by B^^ = Bi + Dne^. This 
would yield /mca = + ^i^o + ^2^0 with ki = Bi{l + 2C12/C11) and /c2 = -0(1 + 2C12/C11) 
instead of eqs. (6,7). Inserting the value of 2C12/C11 for Fe it becomes obvious that D and 
Dn even have a different sign; the sign of D being opposite to the one of the experimentally 
determined Du according to eq. (8). Guo et al. [0] have performed an ab-initio calculation 
of the magnetoelastic properties of epitaxial Co and Ni ffims, and they indeed proceeded 
on this line, i.e. they neglected the second-order term in eq. (1) and instead replaced the 
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constant Bi by a strain-dependent term Si(eo). It should be cautioned that the coefficient 
D obtained from the hnearization of -Bi(eo) is not identical to the second order coefficient 
Dii of the bulk material but it may deviate strongly 
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FIGURES 

FIG. 1. The convergency test for /mca with respect to N^^. The horizontal hnes represent 
the estimated accuracy ±1 ;ueV/unit cell; +, modified tetrahedron method; x, o and □, Gaussian 
smearing with F = 10.2 eV, 6.8 eV and 1.02 eV. The proper convergency using the Gaussian 
smearing is obtained for the smearing parameter where T is from the interval between 6.8 eV 
and 10.1 eV. 

FIG. 2. The calculated magnetocrystalline anisotropy energy density /mca as function of the 
prestrain eo. The solid lines represent fits to quadratic polynomials according to eq. (4). Note that 
ko is indeed negligible; □, LSDA; x, GGA. 
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TABLES 

TABLE L The calculated results in comparison with the experimental data from the thin film 



experiments. 

a[ao] 2g^ Bi[MJ/m3] Dn[GJ/m^] 

LSDA 5.20 1.08 -9.26 1.5 

GGA 5.34 1.13 -2.39 1.1 

exp. 11,1 5.42 1.17 -3 1.0 

exp. II -3.2 1.1 
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